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ABSTRACT 


Astronauts and equipment orbiting the earth 
within the Van Allen Radiation Belts are subjected 
to bremsstrahlung radiation; this radiation is pro- 
duced when electrons are decelerated by interaction 
with matter. The purpose of this paper is to develop 
a capability to evaluate the bremsstrahlung hazard to 
Apollo astronauts; computer programs have been written 
for predicting this component of the total dose. 

These theoretical calculations will be compared in a 
later paper with bremsstrahlung measurements obtained 
from GT-10 . 



BREMSSTRAHLUNG DOSE CALCULATIONS 


FOR APOLLO MANNED SPACECRAFT 

By Joseph W. Snyder, Charles Sterling Portwood, III, 
Donald A. Witt, Timothy T. White, and Alva C. Hardy 


SUMMARY 


Bremsstrahlung, X-rays covering a continuous spectrum, is produced 
during deceleration of electrons as they collide With atomic nuclei. 
Astronauts encounter bremsstrahlung in the Van Allen Radiation Belts, 
where electrons are decelerated in spacecraft shielding; the subsequent 
bremsstrahlung component of total dose is important for mission planning 
and crew safety. The purpose of this study is to develop an ability to 
calculate the bremsstrahlung doses within the Apollo Command and Service 
Modules (CSM) and Lunar Module (LM) during future missions. Specifi- 
cally, the effort has been directed toward the Apollo program, but the 
techniques are applicable to all missions within the Van Allen Belts. 

The computer program, SPW20, was developed to generate dose versus thick- 
ness curves for monoenergetic electrons. The output of SPW20 and that 
of the North American Aviation Geometry Code are used in conjunction with 
MSC-Radiation and Fields Branch dose calculational programs and the MSC 
Orbital Dose Code to compute specific mission doses. 


THEORY 


Introduction 

Bremsstrahlung is energy lost in the form of X-rays by charged par- 
ticles as a result of inelastic collisions with atomic nuclei. The 
energy loss is significant for very light particles such as electrons. 
The total energy of bremsstrahlung emitted by a particle of initial 
energy E q which is completely stopped in an absorber of atomic num- 
ber Z is given by 

I c = kZE o 2 (1) 

where k is taken to be 0.0004 /MeV (ref. 1). 
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The electron range is determined primarily by ionization energy 
losses since bremsstrahlung accounts for only a very small fraction of 
the electron energy. The range of an electron is approximately propor- 
tional to its energy; thus, the electron slows down nearly uniformly 
and most of the bremsstrahlung is produced near the surface where the 
electron enters. 


First Approximation 

The following assumptions are made to simplify calculations: 

p 

1. The electron flux F (electrons /cm /sec) is omnidirectional, 
isotropic, and monoenergetic with energy E o (MeV). 

2. The shield is a spherical shell, outer radius p, thickness T, 
with a dose point at the center. 

3. An electron hitting the shell gives off I Q (MeV) in brems- 
strahlung at the surface, 

4. The bremsstrahlung is created isotropically. 
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The first objective is to calculate the energy flux S (MeV/cm“/see) 
at the dose point, which can easily be converted to dose. The current 
of electrons hitting the shell is defined as 

J s = /j(ft)A eff (ft)dft (2) 


where j is the flux' in a given direction, & e ff ihe cross-sectional 

area in that direction, and ft is solid angle. Under these assumptions, 
this 'becomes 




dft » Ttp F(elec/sec) 


(3) 


Since I is the energy from each electron, the energy flux S' 


of the spherical X-ray source is just 


J I 


S’ = 


8-0 - i pi (, 


l+TTp^ 


A 


MeV/ cm /sec 


■) 


(4) 
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If the average attenuation coefficient over all bremsstrahlung energies 
for the shield is y, then 

S * S'e" MT ■ jj- I , I o e" liT (5) 

Depth Production 

The first modification to the simple model above was to account for 
bremsstrahlung being produced at various depths in the shield, hence 
being attenuated less. Since in a spacecraft configuration p > > T 
and p > > R, the electron range, any small section of the shield may 

be treated as a plane slab. The fol- 
lowing notation will be used (see 
fig. 1): 0 = angle of incoming elec- 

tron with the normal to the shield, 
r * distance along the electron path, 
x * distance along normal to shield, 
x = r = 0 at the surface. It will 
be assumed that T > R (although this 
only changes the limit on some inte- 
grations) and that the electron trav- 
els in a straight line stopping at 
r » R. The next objective will be to 
calculate the profile dl/dx which 
describes the bremsstrahlung produced 
per unit depth. 

By assuming the electron range 
proportional to its energy, the en- 
ergy of an electron is 

Figure 1.- Depth production 

geometry, E = jyi . r/R ) (6) 

and the bremsstrahlung it will give off is 

I = kZE 2 (T) 

Then, since r = x sec 0 (see fig. l) 

S< 9) * (i) (- §) fc) - < 2kZE) (t) (sec 6) 
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|j(e) • -g2(l - r/H)(te« 8) 


Now consider some small area A (cm 4 ) on the surface. The current 
of electrons hitting A is 


m a ![• /*2tt /2 »p 

J a * cos 6 dfl * ^ / d <J> / sin 6 d(sin 6) * -jj- (9) 


Averaging over all particles 


^-f^(e)j(n)A ef . f (a)d « 2 J .^^-(e)cos e sin e]< 
*(i\ 

/ cos l R 1/ v 

' '^1 - ^ sec 0 y sec 6 cos 6 sin 0 d0 


4l o /*cos 


Ul 0 r 008 ' 1 ( R )r 

= R Jo l 


sec 0 I sec 0 cos 0 sin 0 d0 


‘sin 6 d 8 ♦ | d(cos „ e) 
R cos 0 


hr 

dl = o 
dx R 


One may readily show that I 


/* R /dl\ 

J o as 


( 10 ) 


as it should. Now, instead 


of one spherical source, there are many infinitesimal ones having an en- 
ergy flux 


dS' = 


j s di 


Mi )* 1 


4tt(p - x)‘ 


since p > > T > x 


( 11 ) 


Finally, the energy flux at the dose point is 

■ / e ‘ U(T ' X)dS ' * i | e dx (l - | + | in I | |)dx (12) 


which was used in our actual, computations. 
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Anisotropic Correction 

The assumption that hremsstrahlung is produced isotropically is 
not very good for relativistic electrons (E > 1 MeV). The angular dis- 
tribution (ref. 2) is given by 


dl 

dft 


K sin 2 a 



(1-3 cos a)° 


where K is a proportionality constant, o is the angle from the elec- 
tron direction and 3 is the electron velocity divided by the speed of 
light. The first step chosen to correct for this anisotropic distribu- 
tion was to define a correction factor which is the ratio of the fraction 
of bremsstrahlung produced in the inward direction (0 < 0 < tt/2) in the 
actual case to that in the isotropic case. In the actual case we can 
write the fraction inward for any given 0 as 



In evaluating this integral, it is convenient to use an a-i|/ cy- 
lindrical system where a = 0 is the electron direction and ip is the 
azimuthal angle around that direction. The meanings of x and 0 are 
the same as before. The Z-axis is in the plane defined by the X-axis 
and the electron direction, with X-axis lying at ip - 0. The Y-Z plane 
is parallel to the shield surface. The integral can be broken into 
three parts: 

a. If 0 < a < tt/2 - 0, any angle ip is pointed inward. 

b. If tt/2 + 0 < a < tt, any angle ip is pointed outward. 

c. If tt/2 - 0 < a < tt/2 + 0, only the range -tj> m < ip < to is 

pointed inward, where ^ is the value of ip where the X-ray ( y) is 

, emitted parallel to the surface (see fig. 2), 

From figure 2 

sec t|j = -sec fir-ip \=- — = -— *— = - tan a tan 0 (l4) 

ni \ a ha 


o 



Writing equation (13) according to these parts with part (b) = 0 


/•** r ( ' /2, - e di . . 4 /•*»., r ( ’ /2)+e di . . 

Jo d * Jo sin 0 d “ * J-* d * J(i./2)-0 da sin 0 d “ 


r 2 \, r di , . 

jo d * Jo 3n sin 0 d0 


(w/ 2)-8 * r(ir/ 2)+6 

~ sin a do + - /, - sin a da 

dO tt «/(n/ 2 )-e dfl 

f dl < * 

Jo dn sin 0 d “ 


(15) 



In the isotropic case the inward 
fraction is obviously 1 / 2 , so the 
correction factor can then be writ- 
ten as 


f = 
c 


/**/ 2 


u r " /2 

*-J f(e)de 


(16) 


0 A . . , . . , which was used in the actual compu- 

Figure 2.- Anisotropic bremsstrah- ... ^ 

, ... . tat ions, 

lung production geometry. 


Dose Conversion 


Having the energy flux S from equation (32), one may convert to 
dose rate by multiplying by o, the cross section of tissue for the 
bremsstrahlung (averaged over energy)* and an appropriate unit conver- 


2 .ft 

sion factor. With o in cm /g, the conversion factor is 1.6 * 10” , 


so that the dose rate is 


77 = 1.6 x 10 ^of S( rad/ sec) 
at c 


(IT) 


lii 
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SPACECRAFT SHIELDING BREAKDOWN 


The geometry code, developed by North American Aviation, Inc., used 
by MSC in radiation-dose computations utilizes a ray tracing method to 
calculate shield thickness as a function of solid angle for any point 
within a spacecraft or other shielding configuration. The spacecraft 
geometry is defined with all components represented by hexahedrons , cyl- 
inders, spheres, hemispheres, cones, truncated cones, ellipsoids, and 
truncated ellipsoids. These elemental volumes are defined as being 
either solid or void shields with each having its material composition 
and density identified. Complex spacecraft components can be described 
by combinations of different elemental shapes , with void shields deleting 
the interior or other portions of the solid shields. Figure 3 illus- 
trates the description of one such component. The computer program 
calculates thicknesses of the intersected shields along a given ray, 

2 

converts to equivalent thickness in gm/cm of a predefined standard ma- 
terial, and sums the results to obtain the total shield thickness for 
the ray or solid angle. 

The program -considers self-absorption of the crew by the inclusion 
of up to six phantom men similar to the one shown in figure 4. These 
phantom men are used optionally, and may be placed in the spacecraft in 
various standing or sitting positions. 

The program was designed primarily for use in computations of the 
attenuation of high-energy protons and other heavy particles; however, 
it has been adapted for use with electrons and other types of penetra- 
ting radiation. 


BREMSSTRAHLUNG DOSE CALCULATIONS 


In order to support the manned spaceflight program with information 
concerning the radiation hazards of a particular spacecraft trajectory, 
the MSC Orbital Dose Code (MODC) was developed to utilize existing 
environment and dose-calculational programs in obtaining estimates of 
the dose received while in earth orbit. The radiation hazard which 
affects earth-orbital manned spaceflight is due in part to naturally 
occurring electrons and protons which are trapped by the earth’s mag- 
netic field. The problem of trapped radiation was intensified by the 
Starfish high-altitude nuclear detonation of July 9» 1962, which arti- 
ficially injected electrons into the trapping magnetic field. The 
intensity of this artificial radiation is expected to decrease exponen- 
tially with time, but is still of significance when assessing radiation 
hazards. In addition to the radiation from primary electrons and pro- 
tons, bremsstrahlung radiation is also considered. 



Description of spacecraft component by combinations of elemental shapes. 
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The bremsstrahlung dose for a given spacecraft and trajectory is 
determined in the MODC by calculating the instantaneous dose rate for 
numerous points along the trajectory and then integrating numerically 
with resp set to time. The time interval between points must be small 
enough so that the integral can be adequately computed by the trapezoi 
dal rule. The calculational method used in determining bremsstrahlung 
dose is given by 


r E 

I max 

Bremmsstrahlung dose = 57 F^Dk^Dt^ I (dN/dE).jD(E)dE 

1 1 1 1J E 

co 


where 

F^ = the integral electron flux associated with spatial 
location of the ith point on the trajectory 

Dk^ = the fraction the electron flux has decreased due to 
an exponential decay with time of the artificially 
injected electrons for the ith point 

Dt^ = the time interval between the ith and the ith-1 point 

E = upper limit of integration 

max / * 

E co = the average minimum energy of the incident electrons 
sufficient to reach the dose point 

(dN/dE). = the differential of the normalized integral electron 
1 energy spectrum for the ith point 

D(E) = bremsstrahlung *dose as a function of incident electron 
energy for a given spacecraft geometry 

dE = differential electron energy 

The electron flux and energy spectra used by the MODC are taken from 
a recent trapped radiation environment model (ref. 3). The Dk^ s are 

calculated' using the MSC "Decay Program" and the D(E) array was generated 

by MSC - Radiation and Fields Branch dose-calculi tional programs,^" which 
: utilize the results of SPW20. 


1 See section SPACECRAFT SHIELDING BREAKDOWN. 
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CONCLUSIONS 


The computer program SPW20 (see appendix A) produced dose versus 
thickness curves for monoenergetic electrons. The results of this pro- 
gram are shown in figure 5, normalized to one electron per square centi- 
meter. The incident electron energy is at the right of each curve. 

To compare our results with Keller and Pruett (ref. 4), the program 
BREMDOSE (see appendix B) was written. The results of this program are 
shown in figure 6 and are within a factor of 1.5 of similar graphs in 
reference 2. The soft spectra is given by 


F(E) = 1.724 exp 



which is somewhat similar to Van Allen Belt spectra. The fission spec- 
trum is given by 


F(E) = 0.71 exp (-0. 5752 - 0.055E 2 


Because of good agreement between our results and others, using the 
graphs in figure 5 to estimate the bremsstrahlung dose for Apollo seems 
justified. 
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